******************
**** Figure 9 ****
******************

* Run table b15 to get estimates R *
do "$dofiles/table b15" 


coefplot R, keep(treatmentalldummy2 treatmentalldummy3 treatmentalldummy4 treatmentalldummy5 treatmentalldummy6 treatmentalldummy7) vertical yline(0, lcolor(black)) recast(bar) barwidth(0.9) lcolor(gs0) fcolor(*.5) color(gs6) ///
ciopts(recast(rcap)) citop  coeflabel(treatmentalldummy2="P-Bas" treatmentalldummy3="N-Bas" treatmentalldummy4="P-Req" treatmentalldummy5="P-Attr" treatmentalldummy6="N-Req" treatmentalldummy7="N-Attr") format(%9.2f) ///
msymbol(t) levels(95) graphregion(color(white)) addplot(scatter @b @at, ms(i) mlabel(@b) mlabpos(2) mlabcolor(black)) ///
ylabel(-0.15 "-15%" -0.1 "-10%" -0.05 "-5%" 0 "0%" 0.05 "5%" 0.1 "10%" 0.15 "15%" 0.2 " ", angle(0)) ///
ytitle(ITT Effect on Collection Efficiency) title(Overall Treatment Effects, color(black)) subtitle("June 2017 - May 2018") ///
groups(treatmentalldummy2 & treatmentalldummy3  = "{bf:BASIC}" treatmentalldummy4 & treatmentalldummy5 = "{bf:PIDA}" treatmentalldummy6 & treatmentalldummy7 = "{bf:NIDA}")



graph export "$filepath/output_logs/figure 9.pdf", replace
